Infinite matrix product states, Conformal Field Theory 
and the Haldane-Shastry model 
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We generalize the Matrix Product States method using the chiral vertex operators of Conformal 
Field Theory and apply it to study the ground states of the XXZ spin chain, the Ji — J2 model 
and random Heisenberg models. We compute the overlap with the exact wave functions, spin-spin 
correlators and the Renyi entropy, showing that critical systems can be described by this method. 
For rotational invariant ansatzs we construct an inhomogenous extension of the Haldane-Shastry 
model with long range exchange interactions. 



In recent years the study of one-dimensional quantum 
lattice systems has received an enormous impetus due 
the invention of novel numerical algorithms, like DMRG 
[if and its extensions (for a review see Q). The key for 
the success of those methods is the careful treatment of 
the complex structure of the low energy states, which is 
based on the behavior of the entanglement entropy at 
zero temperature Q. For a gapped system, the entan- 
glement entropy, Sl, of a subsystem of length L con- 
verges towards a constant value, independent of L jij. 
For critical systems, however, Sl grows as the logarithm 
of L, were the proportionality factor is related to the 
central charge c of the underlying Conformal Field The- 
ory (CFT) 5]. This behavior of S L implies @] that the 
ground state of the system is well approximated by a 
matrix product state (MPS) Q, a state characterized in 
terms of certain matrices A s . This explains the enor- 
mous success of DMRG, since it can be understood as 
a variational method with respect to MPS of bounded 
dimension ||. However, the unbounded increase of the 
entanglement entropy makes DMRG much less accurate 
for critical systems since the size of the matrices char- 
acterizing the MPS must inevitably grow with the size 
of the system. Thus, in order to properly describe the 
low energy states of such systems MPS with infinitely 
dimensional matrices A s are required. 

The purpose of this letter is to construct such infi- 
nite dimensional MPS (IMPS for short) and their ap- 
plication to ID spin chains. We will replace the fi- 
nite dimensional matrices A s by chiral vertex operators 
of a CFT. These operators act on infinite dimensional 
Hilbert spaces, which acquire the meaning of ancillary 
space where entanglement and correlations are trans- 
ported. The huge enlargement of the ancillary space will 
allow us to describe both critical and non-critical ID sys- 
tems on equal footing. In particular, we will consider the 
XXZ and J1-J2 models to illustrate the accuracy of our 
method. 

Whereas standard MPS have, in general, a complicated 
form if written in the spin basis, the IMPS that we shall 
consider below have a rather simple form given by Jas- 
trow type wave functions. There exists many important 



examples of that type wave functions, such as the Laugh- 
lin of the FQHE 9], the Haldane-Shastry of the inverse 
square Heisenberg model [lOj and the related Calogero- 
Sutherland state of hard core bosons \u\ . A particular 
example of IMPS, which has spin rotational invariance, 
will coincide with the Haldane-Shastry (HS) wave func- 
tion. This case is particularly interesting because there 
is an enhanced symmetry described by the SU(2) Wess- 
Zumino-Witten model at level 1. Using the Ward iden- 
tities of this CFT we shall find an integrable extension 
of the HS Hamiltonian parameterized by a large number 
of coefficients. We will then use this extended HS model 
to study the entanglement entropy of spin chains with 
frustration and random couplings, finding a very good 
agreement with a law recently proposed in 2^, 21 1. 

Let us consider a spin-1/2 chain with N sites, and local 
spin basis \si), where Sj = ±1 (i = 1, . . . , N). A generic 
state in the Hilbert space of this spin chain is given by 
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A MPS state is an ansatz of the form 



i>( Sl ,...,s N ) = (u\A<g...AW\v) 



(1) 



(2) 



where Ag) and (it|, \v) are D-dimensional matrices and 
row and column vectors, respectively. The number or 
parameters needed to specify a MPS scales as 2D 2 N, 
which is much smaller than the number of components 
of a generic state (JlJ, namely 2^. To define a IMPS 
we shall replace the matrices As} by the chiral vertex 
operators of a CFT, thus effectively working with D = 
00. Let us choose for simplicity a gaussian CFT with 
central charge c = 1 . Chiral vertex operators are normal 
ordered exponentials of the chiral bosonic field ip(z) [13], 
which we can be used to define 



A^ = 



Xs 



(* = !,..., AO 



(3) 



where z\ are complex numbers, a a positive real number 
and Xsi = L Si for i even and odd respectively For 
the vectors (u\, \ v), we take the outgoing and incoming 
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FIG. 1: (a) Overlap of the variational ansatz ([5]) with the 
exact GS wave function of the Hamiltonian for N = 20 
spins; (b) Optimal value of a. The dashed line shows the 
curve A = — cos(27ro). 



vacuum states of the gaussian CFT. Equation J2]) be- 
comes then the vacuum expectation value of a product 
of vertex operators which is given by 



JY 



V>(S1, . 



, s N 



(4) 



where 8 S = 1 if X)j=i s i — ^ an d zero otherwise. This con- 
dition can be traced back to the conservation of the U(l) 
charge of the gaussian CFT. Equation® associates the 
charge of the vertex operator a£] to the local spin s, , so 
that the total third component of the spin, S z = J2i s i/2, 
vanishes. We will use the IMPS variationally, with Zi 
and a as variational parameters. The wave function ^ 
scales with an overall factor under a Moebius transfor- 
mation Zi — > (azi + b)/(czi + d) Vi, so that the ansatz 
only depends on N — 3 z-parameters. The latter param- 
eters must not be identified with the position of the local 
spins Si. An exception is the wave function associated 
to translationally invariant spin chains where we shall 



choose Zf, 
have 
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where Sat + i = Si for periodic boundary conditions (we 
take N to be even). This model displays three differ- 
ent phases at zero temperature: i) a gapped AF phase 
(A > 1), ii) a critical phase (—1 < A < 1) and iii) a ferro- 
magnetic phase A < — 1. To describe the GS of the XXZ 
chain we shall use the translational invariant state (0, 
which favours antiparallel nearest-neighbour spins. We 
expect this ansatz to adequately describe the A > — 1 
regimes. The only variational parameter left in |5]) is 
a, which we fix by numerically minimizing the energy. 
In Fig. [1] we plot this overlap as a function of A, for 
N = 20. As anticipated above, the overlap is rather 
poor in the ferromagnetic regime but surprisingly good 
for A > — 1, even above 99%. Moreover, the overlap is 
exactly one at A = and —1, where a = and 1/4 re- 
spectively. This indicates that the corresponding IMPS 
states coincide with the exact GS. To prove this result, 
one applies the unitary transformation U = IIrodd(^f ) 
to ([7|) which flips the signs of the XY exchange inter- 
actions. For A = — 1, one obtains the isotropic fer- 
romagnetic Heisenberg Hamiltonian, whose GS is fully 
polarized with S z — N/2. Applying the lowering op- 
erator (S~) N / 2 to this state, and undoing the unitary 
transformation, one recovers (0 for a = 0. In the case 
A = 0, the transformed Hamiltonian describes free hard- 
core bosons, whose GS is the absolute value of a Slater 
determinant which yields 
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For later purposes, it is convenient to express this spin 
wave function using hard-core boson variables. Mapping 
the spin up (down) states into empty (occupied) hard- 
core boson states, one finds 

iP( ni ,...,n N/2 )Ke^ n * J] ( S m nn \ T nj) ) (6) 



where n\, . . . , n N / 2 are the positions of the bosons on the 
lattice. Undoing the U transformation, one recovers the 
state ([6]), with a = 1/4. For the isotropic XXX model, 
A = 1, the maximal overlap is obtained for a = 1/2. 
The corresponding wave function (0 does not yield the 
exact GS, but its hard-core version © coincides with the 
GS of the Haldane-Shastry model, which belongs to the 
same universality class that the XXX chain (see below 
for a more detailed discussion). 

Collecting the previous results we see that the IMP 
ansatz ([5]), in the range < a < 1/2, corresponds to 
the critical XXZ chain. To confirm the critical proper- 
ties of this ansatz we have computed the Renyi entropy 
— — logTr p 2 L , where is the density matrix of ^ 
restricted to a subsystem of size L. The computation is 
performed as follows. We rewrite 



where = 1, . . . , N denote the positions of N/2 hard- 
core bosons in the lattice. 

1.- The XXZ spin chain: The Hamiltonian with 
anisotropic coupling A is given by 



= \{n,m\il})\ 2 \(n ,m'\ 



n.n' ,m,m 



Hxxz - ^ ( s i s i+i + s f s f+i + &S*Sf +1 ) (7) 



, {ip\n', m)(ip\n, ml) 
(ip\n, m) {ip\n' , m') 

(9) 

where \n) (and \n')) is an orthonormal basis in the space 
of the L spins, and \m) (and \m')) another corresponding 
to the rest of spins. This expression can be easily evalu- 
ated by using two independent spin chains and perform- 
ing the additions using MonteCarlo techniques. In Fig. 
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FIG. 2: (a,b,c) Renyi entropy sf ] for a chain with N = 200 
sites and several values of a £ (0, 1/2]. The dashed line shows 
the CFT prediction for c — 1. (d,e) Correlator C" for several 
values of a. In (d) the stars indicate the exact result. 



EJa,b,c) we plot exp[4(Sf } - S%} 2 )/c] for 2V = 200 and 
several values of cv. For a CFT with central charge c = 1 
one expects this quantity to be sin(7rL/7V) which is 
also plotted (dashed line). The numerical results agree 
(in average) very well with this prediction, confirming the 
criticality of the IMPS for a < 1/2. Note that the os- 
cillations are not a feature of the numerical calculations 
but they are intrinsic in the model. For a > 1/2, we have 
checked that the Renyi entropy saturates to a constant 
value independent of L, in accordance with the gapped 
character of the XXZ spin chain with A > 1. 

Another quantities that can be easily computed using 
MonteCarlo techniques are spin-spin correlators C° a = 
( a n (J o) ( a = x iU: z )- FigHKd,e) shows the correlator 
C* z for a = 1/2,3/4. In the case a = 3/4, the cor- 
relator exhibits anti-ferromagnetic long range order, as 
expected for the gapped A > 1 regime. For a — 1/2 we 
can compare with the exact result (stars in the figure) 
C* z — (— l) n Si(7rn)/(7rn) where Si(x) is the sine inte- 
gral function. For this value of a, the system is isotropic 
(see below), so C° a is independent of a, and its expres- 
sion was first obtained from a Gutwiller proyection of 
the one-band Fermi state 14[. This result motivated 



..(2) 



Haldane and Shastry to construct their inverse square ex- 
change Hamiltonian 10]. These authors also noticed that 
C% x coincides with the one-body density matrix of the 
Calogero-Sutherland model of a gas of hard core bosons 
[ill ]. The GS of the latter model is a continuum version 
of © for a = 1/2. For other values of a E (0,1/2], 
the asymptotic behaviour of the one-body density ma- 



trix correlator is Ci, 



I15H , which can be com- 



pared with the exact scaling in the critical XXZ model 
C%* ~ n^ 71 (A = — cos7r?7) [3]. The comparison of 
these wo results yields, A = — cos(27ra), which correctly 
reproduces the cases: a = 0, 1/4, 1/2 (see Fig. 0Jb)). 



2.- Ji — J2 Heisenberg model: We use the IMPS ansatz 
for this frustated antiferromagnet model described by the 
Hamiltonian 



H 



N 



(JlSj • S,+ l + J2&i ■ Si+2) 



(10) 




FIG. 3: Ji — J2 model (a) Values of 0o and 5 for which the 
energy is minimal; (b) Overlap of the IMPS state with the 
exact GS. The values of the z n in the complex plane are qual- 
itatively depicted for each region. 

The phases of this model are: i) critical c = 1 phase 
(0 < J2 < J2, c ~ 0.241, ii) spontaneously dimerized 
phase ( J2,c < J2 < Jmg = 0.5) and iii) incommesurate 
spiral phase (J2 > Jmg) [13] ■ ^mg is the Majumdar- 
Gosh point whose two GS are the dimer configurations. 
Since this model is isotropic we shall choose the ansatz 
with a = 1/2 (see [HI for a related anstaz). For the 
choice of the Zi parameters we distinguish between even 
and odd sites 



Z2n 



o &-i$0 + —. 



Z2n+1 



-5+i0o + i 



1 N 

, n = 1, . . . , — 
' 2 



where 4>q € [0,7r/iV] and 5 are variational parameters, 
which for a translational invariant state are ((j>o,8) = 
(7r/iV, 0). These parameters are found minimizing the en- 
ergy, and their value (a) as well as the overlap with the 
exact GS wave function (b) are shown in Fig|3] as func- 
tions of J2I J\- At J2 ~ J2,c, the parameter <fio departures 
from tt/N, which reflects the spontaneous dimerization 
of the system. At the MG point, <f>o = $ = and the 
IMPS is a linear combination of the dimer states. For 
J2 > Jmg the parameter 5 departures from zero. Fi- 
nally, at J2/J1 ~ 1, one finds S = tt/N, while S increases 
stcadly with an overlap approaching one, meaning that 
the two chains become increasingly decoupled. A picto- 
rial representation of the variational parameters is also 
shown in figdjb). 
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3.- Inhomogeneous Haldane-Shastry Model: Each (fi- 
nite) MPS is the GS of a finite range parent (frustration 
free) Hamiltonian Hence one may expect the IMPS 
parent Hamiltonians to have a range comparable to the 
size of the system. We will show that this is the case 
for the case a — 1/2 and arbitrary z n . Remarkably, the 
corresponding parent Hamiltonian has two-spin interac- 
tions only and is a generalization of the HS model. First 
notice that the wave function ((4j) is the chiral confor- 
mal block of N primary fields, of spin 1/2 and conformal 
weight h = 1/4, in the SU{2) 1 WZW model. These 
primary fields are nothing but the chiral vertex opera- 
tors ([3]) for a = 1/2. The chiral conformal blocks of 
the WZW model, on a cylinder, satisfy the Knizhnik- 
Zamolodchikov equations, 



+ 2 d_ 



</>cyl(z)=0 (11) 



where i/) C yi(z) = \\f z i^ip( z ) is a conformal transforma- 
tion of the wave function |4]) from the complex plane into 
the cylinder. Now, using eqs pip and (H]) one can prove 
that ^cy^z) is an eigenstate of the Hamiltonian 



with eigenenergy 



WjjjCj - Cj) 

12 



N(N+l) 
16 



Si • S j 



where 



(Zi + Zj)/(Zi - Zj),Ci = £ 



(12) 



(13) 



Taking 



z n = e 2TrTl / Ar , the parameters c, vanish and (fT2")) becomes 
the HS Hamiltonian whose GS is indeed (0. For other 
choices of z n , we have checked numerically that ipcy\ is the 
GS of (|12p . The Hamiltonian ([12")) is an inhomogeneous 
generalization of the HS model, which seems to be inte- 
grable due to the existence of higher conserved quantities 
that can be derived from the KZ equation pT|) . The uni- 
form HS model has a huge degeneracy in the spectrum 
that can be explained by a Yangian symmetry. This de- 
generacy is broken in the non uniform case. 

We can use the above construction to study the en- 
tanglement entropy of random models. The scaling of 
the von Neumann entropy for the random AFH model is 
| log 2 log L [HI , and it is conjectured that the same law 
holds for the Renyi entropy [SEl]]. We take the IMPS 
ansatz g]) with a = 1/2 and z n = e 2 "( n +^)/ N , and ran- 
dom choices of <j) n uniformly distributed in [— 5/2, 5/2]. 
This corresponds to having the inhomogeneous HS model 

(2) 

with random couplings. We have computed averaged 
over realizations and plotted it in Fig. [4] for N = 1000. 
We have subtracted | log 2 log L to compare with the pre- 
diction 2(3, 21 1 which should yield a horizontal line at 0. 



For 5 = 0.75, 0.5, if we ignore the oscillations, we find 
a very good agreement. Even though for 5 = 0.1 the 
result seems to fit better with the formula |logL, the 
inset shows that for long distances the previous result 
seems to apply. Thus, the larger the 5 the more valid the 
prediction is for shorter blocks. 




(21 

FIG. 4: Renyi entropy S L as a function of it = 
log[sin(7r_L/A r )]. We have substracted from the ent ropy 
log(2)xL/3 in order to compare with the prediction [20L [211] . 
Upper (red), medium(green) , and lower (blue) curves corre- 
spond to 5 = 0.75,0.5,0.1. The dashed line corresponds to 
xl/4:. In (b) we have averaged over consecutive points and 
plotted for low values of xl- 



In this letter we have proposed an infinite dimensional 
version of the MPS using the chiral vertex operators of a 
c = 1 CFT. This generalization allow us to study the en- 
tanglement properties and correlators of critical and non 
critical spin chains. For isotropic spin chains the CFT is 
the SU(2)i WZW model and this fact allow us to con- 
struct their associated parent Hamiltonians. They are 
given by an inhomogenous generalization of the Haldane- 
Shastry Hamiltonian, that seems to be integrable. Our 
construction can be generalized in principle to any CFT. 
For rational CFTs, such as the WZW models, we ex- 
pect that the conformal blocks may provide degenerate 
GS with non-abelian statistics encoded in their braiding 
matrices. 
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